Spin-ordering in S = 1 anisotropic Heisenberg models with nondiagonal spin exchange 
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The properties of S = 1 anisotropic Heisenberg models with 
nondiagonal exchange between axial and planar spin compo- 
nents are investigated using Monte Carlo techniques. The 
quantum nature is taken into account in a semi-classical ap- 
proximation. The ordering of the spins when applying an 
external field with axial and planar components is discussed. 
It is argued that the quantum nature of the spins and the 
nondiagonal exchange may explain the peculiar shape of the 
magnetic specific heat of FeBr2 as well as the weakly first- 
order phase transition observed in the same compound when 
a tilted field is applied. 
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I. INTRODUCTION 

Metamagnets of Ising type display a field-induced 
phase transition from an antiferromagnetic phase to a 
(saturated) paramagnetic phase, the transition being of 
first order at low temperatures and of second order at 
high teniperatures.EJp The first and second order transi- 
tion lines meet at a tricritical point. The layered anti- 
ferromagnet FeCl2 is generally considered to be the text- 
book example of this kind of antiferromagnets with a 
strong uniaxial anisotropy.EI 

Recently, there has been a renewed int^est in the 
layered hexagonal antiferromagnet FeBr2Jjil3 In the or- 
dered phase of FeBr2 the spins of the iron ions are aligned 
ferromagnetically in the triangular layers, the layers be- 
ing stacked antiferromagnetically along the c-axis. Ad- 
jacent iron layers are separated by two nonmagnetic bro- 
mide planes. |-| 

Investigations of de Azevedo et al.Q revealed noncrit- 
ical anomalies in the antiferromagnetic phase of FeBr2 
below the transition line to the paramagnetic phase, e.g. 
maxima or shoulders in the temperature derivatives of 
the magnetization at constant axial fields. The ingre- 
dients crucial for the existence of these noncritical spin 
fluctuations are the effectively weak ferromagnetic in- 
tralayer couplings and the large number of interlayer 
nearest neighbors (due |-t|a|-jbhe superexchange mediated 
by the bromide planes) .ErElla 

Surprisingly, measurements of the magnetic specific 
heatQ revealed a sharp peak superposed on a broad shoul- 
der or maximum below the transition to the paramag- 
netic phase, in con-tta^t to Monte Carlo simulations of 
Ising metamagnetajljlj showing only a noncritical shoul- 
der or maximum. The peculiar shape of the specific 



heat of FeBr2 was tentatively interpretedQfl as the sig- 
nature of a new phase transition between two differ- 
ent antiferromagnetic orderings of the z-components of 
the spins. This new transition line was supposed to 
result from the decomposition of the tricritical point 
into a bicritical point and a critical end point, a pos- 
sible scenario emergiiig from the mean field treatment 
of Ising metamagnets.EJ Remarkably, however, such a de- 
composition of the tricritical point has not been observed 
in Monte Carlo simulations of three-dimensional Ising 
antiferromagnets £3l3^t3 Furthermore, there are no indi- 
cations of an antiferromagnetic-antiferromagnetic phase 
transition in the magnetometric measurements in an ax- 
ial magnetic fieldoEiJ 

Recent investigations in an external field with an axial 
and a planar component have shown evidence for jumps 
in thcL. magnetization parallel and perpendicular to the 
field.E3 The authors suggest that this transition-like phe- 
nomenon may involve a nondiagonal exchange between 
the axial and planar components of the spins. 

Hitherto, theoretical progress in understanding the 
properties of FeBr2 has mainly been achieved,|b5t.a|nalyz- 
ing simplified spin ^ Ising antiferromagnets.EjETQ'El How- 
ever, the lowest state of an iron ion in FeBc^is a triplet 
consisting of a lowest doublet and a singlet .113 Therefore, 
in Ref. 9, anomalies of the magnetization and the specific 
heat have also been studied in S = 1 models. 

In the present work, I shall present a Monte Carlo 
study of S = 1 anisotropic Heisenberg models which in- 
clude a nondiagonal exchange term between axial and 
planar spin components. The quantum nature of the S 
— 1 spins will be taken into account in a semi-classical 
approximation. Especially, I shall give a possible expla- 
nation of the peculiar shape of the specific heat of FeBr2 
and discuss the spin-ordering in an axial and in a non- 
axial field. 

The paper is organized in the following way: In Sec- 
tion II the anisotropic S = 1 Heisenberg models are in- 
troduced and the semi-classical approximation used in 
the simulations is discussed. Section HI deals with the 
properties of ferromagnetic and antiferromagnetic mod- 
els on a tetragonal lattice. These simplified models are 
studied in order to investigate the effects of (1) the semi- 
classical approximation, (2) different couplings between 
planar and axial spin components, and (3) an external 
planar field component. In Section IV a realistic model 
for FeBr2, both in an axial and in a tilted field, is dis- 
cussed and the findings are compared to the experiments. 
Some of the results of tii^ Section have already been pub- 
lished in a brief report .liZI A short summary concludes the 
paper. 



II. THE MODELS 



The low-temperature properties of FeBr2 (space group 
_ j^, may be-des^pjbed by the effective S = 1 Heisenberg 
HamiltonianliS 1 



^3d 



Ti — 2_^ Jq/3 < —S^Sp + S^Sp + S^Sp 



■ E ^ (^S)' - E [^-^" + ^-^^] + '^ndiJ.z) (1) 



defined on a three-dimensional hexagonal lattice con- 
sisting of a layered system of triangular lattice planes 
stacked along the c-axis. Two different sets of interac- 
tions Jap have-been derived from inelastic neutron scat- 
tering studies :t3ll3 the ferromagnetic nearest neighbor in- 
teraction in the triangular layers perpendicular to the 
c-axis (corresponding to the z-direction), Ji, is weak- 
ened by antiferromagnetic in-layer cauplings, extending 
either up tii next-nearest neighbors JIj J2, or up to third 
neighbors ,tj J3. Every spin is coupled antiferromagneti- 
cally pith the strength J' to ten spins in the neighboring 
layer.Ej The constant 77 describes the anisotropy in the 
exchange interactions. 

The second term in Equation (1) describes a single-ion 
anisotropy with the easy axis of the spin along the z- 
axis, i.e. D > 0, whereas Hx and Hz are the planar and 
axial components of an applied external field. Finally, 
the last term, TindiJxz), is a nondiagonal intralayer ex- 
change between axial and planar spin components with 
the strength Jxzf^ involving only products of pairs of 
spins. In order to illustrate this interaction, which is 
invariant when applying the symmetry elements of the 
trigonal point group, consider two neighboring in-layer 
sites, called and 1, with the same value of y.pThe non- 
diagonal exchange may be written in the formB2l 



Ti-ndiJxz) — E ^afsiJxz) 
(a/3) 



(2) 



where the sum is over nearest-neighbor in-layer sites 
(a/?)|_.For the pair (a/3) = (01) this exchange takes the 
formEj 



Hqi — —Jxz {SqS^ + SqS[) . 



(3) 



The other terms in Equation (2) arc obtained by applying 
the appropriate symmetry elements transforming the pair 
(01) onto (a/3). 

As this model, defined on a hexagonal lattice, is 
rather complicated, it is preferable to study the effects 
of the semi-classical approximation, of the nondiagonal 
exchange, and of the planar field component first in sim- 
pler models. Therefore, I consider also the following S 
= 1 Heisenberg Hamiltonian defined on a tetragonal lat- 
tice: 



'H — —J y ^ Sjjk iSi^ijk + Sij-^-lkj — J 7 ^ SijkSjj k+i 

ijk ij k 

-D E {s0' - E ^4^0 + nndiJxz) (4) 

ij k ij k 

with S = {S'',Sy,S^) and H = (iJ^,0,i/^). Here i, j, 
k (corresponding to the a, b, c-directions respectively) 
label the lattice sites. The coupling in the layers, J, 
is ferromagnetic, whereas the coupling between adjacent 
layers, J', may be antiferromagnetic (J' < 0) or ferro- 
magnetic (J' > 0). The single-ion anisotropy, D > 0, is 
considered in order to obtain ground states with nonvan- 
ishing z-components of the spins. This may alternatively 
be achieved by considering instead an anisotropy in the 
exchange interactions (yielding three-dimensional XXZ 
models). 

For the tetragonal models, I will consider two different 
couplings between planar and axial spin components: (1) 
a term involving only products of pairs of spins, being in- 
variant under the symmetry operations of the point group 
C^v (this term can thus be considered to be analogous to 
the nondiagonal exchange term, see Equations (l)-(3), 
proposed for the description of FeBr2), 



'^ndK'Jxz) ~ ''xz / , \^ijk'^i+l jk ^ijk^i+ljk 
ijk 

^'-'ijk'-'ij+lk 



_ qV qz 

•^ilk^ij + lk 



(5) 



or (2) a term involving products of squares of spins, being 
invariant under the symmetry operations of the space 
group. 



''■ndi'^xz 



— Jxz 2^ 
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where the sum is over nearest-neighbor in-layer sites. 

I consider in the following S = 1 spins where the quan- 
tum nature is treated in a semi-classical approximation: 
the z-component is discretized and can only take the 
values 1, 0, or —1, whereas the spin length is fixed to 

= V2. 



be 



S 



^S (S + 1) = -\/2. Hence, the spins rotate in 

the xy plane like a classical vector. The xy-components 
provide additional energy contributions as compared to 
the S = 1 Ising model, even at low temperatures. Note 
that in Ref. 9 the spin length was fixed at 1, the spins 
thus having a planar component only if 5^ = 0, yielding 
thermal properties similar to the S = 1 Ising model. 

The ground states of the semi-classical antiferromag- 
netic hexagonal and tetragonal models are readily ob- 
tained provided the nondiagonal exchange constant is not 
too large. In small axial fields {Hx = and H^ < H^ = 
-2/VJ', N = 10 for the hexagonal model and TV = 2 
for the tetragonal model), the antiferromagnetic phase, 
where the spins are aligned ferromagnetically in the lay- 
ers, with an antiferromagnetic arrangement between sub- 
sequent layers, is stable. For larger fields (H^ > H^ and 



Hx = 0) the axial spin components S^ order ferromag- 
netically, whereas the signs of the planar components {S^ 
and S'^) still change from layer to layer. Note that in ab- 
sence of an ordering planar field the ground state is in- 
finitely degenerate, as the angle between the planar spin 
components and, say, the a;-axis is not fixed. Keeping 
Hz constant and applying a planar field, Hx, yields a 
stable phase where the xy-components of the spins are 
ordered in a spin-flop phase, in which the magnetization 
per layer in x-direction has the same value in each layer, 
whereas the j/-components have different signs in adja- 
cent layers but the same absolute value. For stronger 
planar fields the y-components finally vanish, leading to 
a ferromagnetic ordering of the planar spin components. 
The effect, at T = 0, of a nondiagonal coupling be- 
tween spin components can be discussed for the tetrago- 
nal ferromagnetic model. For small absolute values of 
the nondiagonal exchange constant (be it J^z or Jj^, 
see Equations (5) and (6)), the ground state is given by 
a ferromagnetic phase where the spins have both axial 
and planar components. Consider first the two-spin ex- 
change with the coupling constant J^z ■ Increasing J^z 
does not change the energy of the ferromagnetic phase, 
whereas the energy of a second phase, consisting of fer- 
romagnetically coupled layers with a stripped pattern in 
every layer, is decreased, until, for | J^^j > V2J + ^H^, 
the stripped phase has, at T = 0, a lower energy than 
the ferromagnetic phase, thus forming the ground state. 
Note that the appearance of a new ground state at the 
threshold value does not depend on the sign of J^z- If 
the second spin exchange with negative coupling con- 
stants Jxz is considered, the energy of the ferromagnetic 
phase is increased when J^z decreases. For couplings 

-r, d+h^+(i-V2)h^ , , ... 

Jxz < 4 — the ground state is given by 

a ferromagnetic phase where the spins have only a pla- 
nar component. A positive coupling constant Jj^ does 
not yield a new stable phase. The ground state of the 
tetragonal antiferromagnetic model in presence of a non- 
diagonal spin exchange can be discussed in similar terms. 
In the following it is always supposed that the nondiago- 
nal coupling constant does not exceed the threshold value 
which leads to a new ground state. 

The thermal properties of the different models were in- 
vestigated by simulating systems consisting oi L x L x L 
Heisenberg spins, with L ranging from 10 to 30. As the 
equilibration proved to be rather slow, the first 3 x 10^ 
Monte Carlo steps per site were usually discarded. Sim- 
ulations for up to 20 different realizations with different 
random numbers were performed in order to improve the 
statistics. Besides the energy E and the specific heat C, 
the different components of the magnetization per layer, 
mx{i), my{i), and mz{i), and, for the antiferromagnetic 
models, related quantities such as the sublattice magne- 
tizations, Ml — (^ml.,my,ml) and AI2 = ('7^^, "i^, to^), 
referring to odd and even layers, were computed. Here, 
the odd, respectively even layers have the magnetization 
ml = +1, respectively ml = — 1 at T = 0. The applied 



axial field points in the "-I-" -direction, i.e. it tends to 
destabilize the even layers. 

The ground state was always used as starting config- 
uration for the simulations. In absence of a planar field 
component, the state with S'" — was chosen among the 
infinity of degenerate ground states. Note that in this 
case one encounters, after initial relaxation, a metastable 
state, in which the system remains, possibly, for a long 
time, the time depending on the size of the system and 
the temperature. The computed planar spin components 
are then supposed to be very close to their values in the 
thermodynamic limit. 



III. THERMAL PROPERTIES OF THE 
TETRAGONAL MODELS 

In this section, the thermal properties of the tetragonal 
ferromagnetic and antiferromagnetic models (see Equa- 
tion (4)) are discussed. For a vanishing external field, 
the Hamiltonian of the antiferromagnetic model can be 
mapped in the usual way onto that of the ferromag- 
netic model. Therefore, I will in the following present 
for Hx = simulations of the ferromagnetic model, the 
antiferromagnetic model being analyzed for Hx 7^ 0. 

Figure 1 shows the temperature dependent specific 
heat and the magnetization obtained for vanishing nondi- 
agonal spin exchange for the ferromagnetic model, with 
J = J', D = 3J', and Hx = 0. The specific heat has 
a two peak structure, see Figure la: the peak at Tc is 
the critical peak resulting from the disordering of the z- 
components of the spins, whereas the low-temperature 
peak at Txy follows from the disordering of the planar 
spin components. The magnetization data, see Figure 
lb, show that the z-components, in absence of a cou- 
pling between planar and axial spin components, are not 
affected by this disordering. 

This two peak structure of the specific heat should be 
compared to the specific heat of the classical anisotropic 
Heisenberg model showing only one maximum. For van- 
ishing single-ion anisotropy the two peaks of the specific 
heat of the semi-classical model merge to a single peak 
located at the critical temperature of the corresponding 
classical isotropic Heisenberg model. 

It is clear from Figure 1 that, due to the discretization 
of S"^ and the presence of the single-ion anisotropy, the 
disorderings of the axial and planar spin components are 
largely decoupled. The z-components behave like S — 
1 Ising spins and disorder at the critical temperature Tc 
(which increases with increasing value of D, yielding in 
the limit D — > -|-cxd the critical temperature of the 3d 
Ising model), whereas the xy-components form a clas- 
sical three-dimensional plane-rotator disordering at the 
temperaturecj Txy = 2.2J' (setting the Boltzmann con- 
stant equal to one). 

Note that the specific heat does not vanish when T 
approaches 0. This is an artefact of the classical nature 
of the planar spin components yielding limy^o C = i. 



In presence of the nondiagonal two-spin exchange with 
coupHng constant J^^ (Equation (5)), the ordering tem- 
perature Tc of the axial spin components is decreased, 
see Figure 2. This results from the coupling, at temper- 
atures above T^y, of the ordered z-components to the 
disordered planar components, yielding additional fluc- 
tuations which lead to a decrease of the ordering tem- 
perature of the axial spin components. The value of Txy 
is not affected by this nondiagonal spin-exchange, as long 
as \J].J\ is not too large, see below. The spin exchange 
involving products of squares of spins (Equation (6)) also 
leads to a decrease of Tc for negative values of J^^. This 
second coupling is not so effective in destabilizing the fer- 
romagnetically ordered z-components, yielding a smaller 
decrease of Tc. 

If the strength of I J^^ I is increased beyond a threshold 
value (being around 0.65 J' for J = J' and D == 3J'), a 
jump in the axial component of the magnetization, Mz, 
is observed, see Figure 3a. This discontinuous change, 
which is induced by the disordering of the planar com- 
ponents, shows up as a sharp peak in the specific heat 
(Figure 3b). At temperatures slightly above Tj;y some 
orderness of the z-components still persists, i.e. M^ ^ 0. 
Increasing further the temperature, M^ increases until, 
at Tc {— 3.1 J' for the parameters of Figure 3), the dis- 
ordering of the z-components finally takes place, yield- 
ing a second peak in the specific heat. This behaviour 
of Mz may be better understood when the temperature, 
starting from the disordered high temperature phase, is 
decreased. At Tcj the z-components order at the usual 
second order phase transition. The increase of the order- 
ing of the S^ spins at temperatures below Tc increases 
the effect of the nondiagonal exchange, yielding, after an 
initial slow rising, effectively a decrease of M^, until, at 
Txy, the planar components order. Note that in this case 
both Tc and T^y are shifted to lower temperatures, the 
shift being larger for larger values of J\^ . 

Interestingly, no jump in M^ occurs when the second 
spin exchange with coupling constant Jj^ is considered. 
For all values of Jj^, yielding the ferromagnetic ground 
state with both planar and axial spin components, see 
above, M^ changes continuously at T^y- 

Figure 4 shows the specific heat of the antiferromag- 
netic tetragonal model obtained when planar fields with 
different strengths are applied and only diagonal cou- 
plings are considered. Increasing the value of tl^ moves 
Txy to lower temperatures, whereas Tc is not changed. In 
presence of a planar field unusual strong finite-size effects 
are observed when the planar spin components disorder, 
making it necessary to simulate at least systems of length 
L = 20 in order to get a reliable estimate of T^y These 
strong finite-size effects are due to the presence of the 
S"^ = states and do not show up for the correspond- 
ing classical plane-rotator, as I checked. Finally, one 
should notice that, in presence of the nondiagonal cou- 
pling (5) with the strength J^^, the total magnetization 
in z-direction again changes discontinuously for values of 
\j\A larger than the threshold value, which seems not to 



depend on the value of H^ 



IV. THERMAL PROPERTIES OF THE 
HEXAGONAL MODEL 

Most of the features discussed in the previous Section 
for the tetragonal models are also encountered if the r|e^ 
alistic hexagonal model (1) for FeBr2 is considered.EZi 
In the following I will not reiterate the discussion on 
the values of the diagonal couplings of the spins (the 
interested reader is referred to Ref. 9) and choose a 
set of parameters, .based on inelastic neutron scatter- 



ing experimentsli-Hij, for which pronounced anomalies in 
the magnetization data and in the specific heat data are 
observedB: Ji = -16.75J', J2 = 0, J3 = -0.29Ji, and 
7? = 0.78. 

Figure 5 shows the specific heat and the sublattice 
magnetizations obtained for strong single-ion anisotropy 
D and vanishing nondiagonal interactions, with H^ = 
— 18 J' and H^ = 0. The temperature dependent specific 
heat C{T) has a three peak structure: the peak at Tc is 
the critical peak resulting from the disordering of the z- 
components of the spins. Going to lower temperatures, 
one encounters first the anomaly due to the noncritical 
spin fluctuations. This peak does not correspond to a 
sharp phase transition and does not show a significant 
size dependengy, in contrast to the critical peak at higher 
temperatures.lj Both peaks are present when computing 
the specific heat of Ising metamagnets and involve solely 
S^ . The third peak, again, results from the disordering 
of the planar spin components, see Figure 5b. 

The presence of the nondiagonal spin exchange moves 
both Tc and the anomaly to lower temperatures, see Fig- 
ure 6a, whereas T^y only decreases, moderately, for large 
values of J^z- The anomaly in the specific heat ap- 
proaches the xy-peak, its height being reduced, when 
Jxz is increased, and finally merges with the T^y peak. 
Again, a threshold value exists for \Jxz\, at about 18 
J', beyond which the total magnetization in z-direction, 
Mz — [ml +7111) /2j changes discontinuously, see Figure 
6b. The disordering of the planar components at Txy 
is then also clearly discontinuous, as can be seen in the 
jump of the order parameter My'P — {m} — mT) /2, see 
Figure 6c. 

Changing the value of the axial field Hz does not 



change Txy if Jx 



0. For nonvanishing nondiagonal 



couplings, however, the position of the xy-peak is shifted 
to lower temperatures when the axial field strength is in- 
creased. 

Figure 7 shows the influence of the single-ion 
anisotropy D on the specific heat for Jxz — 0. Decreas- 
ing Z?, the Tc peak and the anomaly are moved to lower 
temperatures, whereas Txy is slightly increased. A sim- 
ilar behavior is found for Jxz ^ 0, as shown in Figure 
8 for Jxz = 16. 2J' and D = -8.1J'. For this choice of 
the parameters, the noncritical spin fluctuations appear 



at lower temperatures than the disordering of the planar 
spin components. The resulting specific heat has a pe- 
culiar shape consisting of a broad shoulder, the anomaly, 
and a superposed peak, the T^y peak. This shape is rem- 
iniscent of the magnetic specific heat of FeBr2 and will 
be discussed in more details below. 

Applying an additional planar field component, Hx, 
leads to a spin-flop phase in the xy-components at 
T — 0: the y-components of the layer magnetization 
have opposite signs in the different sublattices but the 
same absolute value, whereas the a;-components are the 
same in every layer. As the axial field tends to stabilize 
the odd or "-I-" layers and to destabilize the even or "— " 
layers, the absolute values of mj and rUy are close but 
not identical at temperatures T > 0. 

Figure 9 shows the different sublattice magnetizations 
for Jxz = 16. 2J', D = -8.1J', and H^ = 0.75H,, i.e. 
a field with axial and planar components is considered. 
At low temperatures, the antiferromagnetic ordering of 
the axial spin components and the spin-flop ordering of 
the planar spin components is clearly seen. When the 
system is heated, a drastic change in the y-component 
of the layer magnetization takes place at T^y well be- 
low Tc- For the chosen parameters, the change in the y- 
components is continuous. For values of | Jxz \ larger than 
the threshold value, e.g. J^z = 21.6 J', this change occurs 
through a jump in the layer magnetization leading to a 
first-order transition. At first it may be surprising that 
the j/-components (and also the x-components) of the 
sublattice magnetization are not equal at temperatures 
above T^y and below Tc, see Figure 9. In fact, this is 
explained by the finite value of the axial field Hz which 
tends to destabilize the z-components in the "— " lay- 
ers, thus leading to different sublattice magnetizations in 
the different layer types. As a result the y-components 
do not vanish at the phase transition at T^y, but take 
a small value which is different for the two layer types. 
Only when the z-components of the sublattice magne- 
tizations are equal, at T^ do the j/-components vanish 
and the x-components take the same value. This effect 
is also observed when applying a field with axial and 
planar components to the tetragonal antiferromagnetic 
model discussed in the previous Section. 

Increasing the value of the planar field moves T^y to 
lower temperatures, as it is also the case for the tetrag- 
onal antiferromagnetic model. Furthermore, the specific 
heat seems to loose its peculiar shape for large values 
of Hx- At least for the temperature resolution used in 
the present study, only one, rather broad, peak is seen 
instead of the shoulder or maximum with a superposed 
sharp peak which is observed for Hx — 0. 

The presented data clearly indicate that the quantum 
nature of the S = 1 spins has to be taken into account 
in a theoretical description of the low temperature be- 
haviour of FeBr2. The semi-classical approach adopted 
in the present work leads to different disordering tem- 
peratures for the planar and the axial spin components. 
Hence, the anomaly and the low-temperature xy peak are 



not intimately related. The modification of, for instance, 
the exchange Jxz or the degree of Ising-like anisotropy 
D changes their respective positions and may lead in 
an axial field to specific heat data having the peculiar 
shape of the magnetic specific heat of FeBr2 , see Figure 
8. One should notice that, in order to obtain this pecu- 
liar shape with the considered set of diagonal coupling 
constants, the single-ion anisotropy has to be decreased 
considerably as compared to the values obtained from in- 
elastic neutron scattering measurements. Nevertheless, 
one must keep in mind that both D and the strengths of 
the diagonal couplings were derived from the experiments 
without itaking a possible nondiagonal spin exchange into 
account £3't2l 

The evolution of the temperature dependent specific 
heat of|-.FeBr2 with increasing axial fields is also of 
interest.!!] For small fields only one peak, the critical peak 
at Tc, has been observed. For larger fields, a shoulder 
with a superposed peak appears at lower temperatures. 
The shoulder evolves into a maximum with increasing 
fields whereas the superposed peak becomes sharper, the 
whole being shifted to lower temperatures. This large- 
field behavior is well rendered in the simulations, as 
shown in Figure 10. Indeed, the shoulder, resulting from 
the noncritical spin fluctuations, changes into a maxi- 
mum when Hz is increased, whereas the xy-peak be- 
comes sharper, the anomaly in the axial spin components 
and the disordering of the planar spin components both 
moving to lower temperatures. Note, however, that in 
the present model the disordering of the xy-components 
also takes place for small fields, see above. 

The magnetization data resulting from the Hamilto- 
nian (1) also compare favorable to the experiments. For 
example, without an applied planar field, the trans- 
verse spin-ordering at Txy does not lead to a discon- 
tinuity of the axial magnetization (if Jxy is not too 
large, see above), in accordance with the magnetometric 
measurementtcl, showing no jump in the axial magnetiza- 
tion despite a sharp peak in the specific heat. An addi- 
tional planar field component yields a low-temperature 
spin-flop phase in the a:j/-ccttHponents, which is com- 
patible with the experiment .IIj The change in the y- 
component in the magnetization per layer may then oc- 
cur through a continuous or a discontinuous phase transi- 
tion, depending on the values of the parameters. Exper- 
imentally, one observes, in a tilted field, a 'weakly first- 
order' transition,ll3 thus being supposedly at the border 
between these two scenarios. 

One should bear in mind when applying the present 
results to FeBr2 that the Hamiltonian (1) may not be 
complete. Indeed, the recent measurements in an nonax- 
ial fielc£3 suggest the presence of in-plane anisotropy in 
FeBr2, resulting, for example, from magnetoelastic cou- 
plings. Thus an anisotropy is not contained in the present 
Hamiltonian. Furthermore, one should recall that the 
treatment of the Heisenberg Hamiltonian is only approx- 
imative, due to the semi-classical approach adopted in 
Section II, leading to a strong decoupling between pla- 



nar and axial spin components. 

Nevertheless, the present study clearly shows the im- 
portance of the nondiagonal spin exchange and of the 
quantum nature of the S = 1 spins in describing the low- 
temperature properties of FeBr2 . Especially, the peculiar 
shape of the specific heat in an axial field may be traced 
back to a phase transition of the planar spin components, 
taking place, rather by chance, close to the anomaly of 
the axial spin components in the antiferromagnetic phase 
of FeBra. 



described if a nondiagonal coupling in the spin compo- 
nents and the quantum nature of the spins are taken into 
account. It seems, in the light of the present study, that 
the anomaly in the axial spin components in the antifer- 
romagnetic phase and the disordering of the planar spin 
components occur in FeBra closely one to the other rather 
by chance, yielding, for instance, the measured peculiar 
shape of the magnetic specific heat. 
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V. SUMMARY 

Recent experimental investigations showed some in- 
triguing transition-like phenomena in the magnetiza- 
tion and in the specific heat of the layered antiferro- 
magnet FeBra. The observed transverse spin-ordering 
in a nonaxial field suggests the existence of a nondi- 
agonal exchaiige between the planar and the axial spin 
componentsJlj 

Motivated by these experiments, the properties of S = 
1 anisotropic Heisenberg models (both on a tetragonal 
and on a hexagonal lattice) including nondiagonal spin 
exchange were investigated by Monte Carlo simulations. 
The quantum nature of the spins were taken into account 
in a semi-classical approximation, where S^ was dis- 

cretized and the spin length fixed to be S = y^S (S -t- 1). 

Applied external fields with axial and planar components 
were considered. 

These simulations showed that, besides the disorder- 
ing of the axial spin components at high temperatures, a 
second phase transition involving the planar spin com- 
ponents takes place at lower temperatures. For the 
tetragonal models, the effects of two different nondiag- 
onal spin exchange terms were discussed. It was shown 
that the nondiagonal coupling involving only products 
of two spins could induce, for large values of the cou- 
pling constant, a discontinuous change of the total axial 
magnetization at the ordering temperature of the planar 
spin components, whereas Mz always changed continu- 
ously when the second exchange term involving products 
of squares of spins was considered. 

For the realistic antiferromagnetic model defined on a 
hexagonal lattice, which is supposed to describe the low- 
temperature properties of FeBra, the jump in the total 
magnetization in axial direction for strong nondiagonal 
couplings was always mirrored by a corresponding jump 
in the planar sublattice magnetizations. Furthermore, it 
was shown that the anomaly of the if^ff-l spin compo- 
nents in the antiferromagnetic phaseB-Eia and the low- 
temperature peak resulting from the disordering of the 
planar spin components are not intimately related. In 
fact, the data suggest that the transition-like features 
in the specific heat (superposed peak on a shoulder or 
maximum) and in the magnetization (jumps in the mag- 
netization observed in a tilted field) of FeBra are well 
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FIG. 2. Specific heat C of the ferromagnetic model as a 
function of temperature for the two different nondiagonal spin 
exchange terms given in Equations (5) (coupling constant J^y) 
and (6) (coupling constant J^y), with J = J' , D = 3,7', and 
Hx = Hz — 0. Open diamonds: J^y = 0.5 J'. Filled squares: 
Jxy = 0.5 J'. The data obtained without a nondiagonal cou- 
pling (open circles) are included for comparison. Systems with 
10^ spins are considered. Only error bars larger than the size 
of the symbols are shown. 



FIG. 1. Monte Carlo data of (a) the specific heat and (b) 
the components of the magnetization as function of the tem- 
perature for the S = 1 ferromagnetic Heisenberg model on a 
tetragonal lattice, with J — J' , D — 3J', and Hx ~ Hz — 0. 
The system size is L = 20. Here and in the following figures, 
the Boltzmann constant is set equal to one. 
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FIG. 3. Monte Carlo data of (a) the components of the 
magnetization and (b) the specific heat as function of the tem- 
perature obtained for the ferromagnetic model on a tetragonal 
lattice, with J = J' , D = 3J', Hx = Hz = 0, and J^z = 0.7J' 
. Systems with 20 x 20 x 20 spins are simulated. 
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FIG. 5. Monte Carlo data of (a) the specific heat and (b) 
the components of the magnetization per layer in odd (full 
symbols) and even (open symbols) planes as function of the 
temperature for the S = 1 anisotropic Heisenberg model on 
the hexagonal model in an axial field Hz = —18 J', with 
D = -32.4 J' and Jxz = 0. Systems with 20 x 20 x 20 spins are 
simulated. Error bars are only shown when they are larger 
than the sizes of the symbols. 
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FIG. 4. The temperature dependent specific heat obtained 
for the antiferromagnetic S = 1 Heisenberg model defined on 
a tetragonal lattice when planar fields with different strengths 
are applied, with J = —J', D — — 3J', and Hz ~ 0. Nondiag- 
onal couplings between axial and planar spin components are 
not considered. The system size is L = 20. Only error bars 
larger than the size of the symbols are shown. 
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FIG. 6. Temperature dependence of (a) the specific fieat 
(the arrows indicate the critical temperatures, Tc), (b) the 
total magnetization in ^-direction, and (c) the order param- 
eter of the planar spin components for different values of the 
nondiagonal spin exchange between the planar and the axial 
spin components for the hexagonal antiferromagnetic model, 
withD = -32 AJ' and H^ = -18J'. Systems with 20 x 20 x 20 
spins are considered. Only error bars larger than the size of 
the symbols are shown. 
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FIG. 9. The components of the magnetization per layer 
in odd (full symbols) and even (open symbols) planes as 
function of temperature in presence of a planar field com- 
ponent, H^ = 0.75//^, with J^^ = 16. 2J', D = -8.1J', and 
Hz = —18 J'. Hexagonal systems with 30 x 30 x 30 spins were 
simulated. 



FIG. 7. The specific heat C as a function of temperature 
for different values of the single-ion anisotropy, with Jxz ~ 
and Hz = —18 J'. The system size of the simulated antifer- 
romagnetic hexagonal model is L = 20. Error bars are only 
shown when they are larger than the sizes of the symbols. 
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FIG. 10. Specific heat C vs temperature as obtained from 
Monte Carlo simulations for different strengths of the axial 
field, with Jxz = 16.2 J', D = -8.1J', and H^ = 0, for hexag- 
onal systems with 20"^ spins. Error bars are only shown when 
they are larger than the sizes of the symbols. 



FIG. 8. Temperature dependent specific heat C of the 
hexagonal model with J^z = 16. 2J', D = — 8.1J', and 
Hz = —18 J', for systems with 20 x 20 x 20 spins. 



